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Abstract. We investigate the possibility and conditions to factorize the entanglement 
evolution of a multiqubit system passing through multi-sided noisy channels. By 
means of a lower bound of concurrence (LBC) as entanglement measure, we derive 
an explicit formula of LBC evolution of the iV-qubit generalized Greenberger-Horne- 
Zeilinger (GGHZ) state under some typical noisy channels, based on which two kinds 
of factorizing conditions for the LBC evolution are presented. In this case, the time- 
dependent LBC can be determined by a product of initial LBC of the system and 
the LBC evolution of a maximally entangled GGHZ state under the same multi-sided 
noisy channels. We analyze the realistic situations where these two kinds of factorizing 
conditions can be satisfied. In addition, we also discuss the dependence of entanglement 
robustness on the number of the qubits and that of the noisy channels. 
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1. Introduction 

Quantum entanglement, as a type of distributed or nonlocal coherence among several 
quantum subsystems, have been more and more recognized as an indispensable resource 
in realizing many intriguing quantum information processing and quantum computation 
p. Realistically, however, the unavoidable coupling of the entangled system to 
environment will lead to a destruction of the necessary entanglement. Therefore, a 
deeper understanding of entanglement dynamics is of great importance not only in 
the foundation of quantum mechanics [2J but also in the rapidly developing quantum 
technologies pp. 

The usual way in studying entanglement evolution is first to deduce the evolved 
state of the system and then to calculate its entanglement [3j HI El El [T_l El [9j [TO] . 
However, as far as multipartite systems (or higher dimensional systems) be concerned, 
it is generally very hard to solve the state equation and therefore the evolution 
of entanglement can be determined only in very special cases. Different from the 
aforementioned state-evolution technique, it was found [UJ for any two-qubit system 
with only one qubit being subjected to a noisy environment, i.e., under one-sided 
noisy channel, the evolution of system's entanglement in terms of concurrence [T_2] 
can be completely determined by the product of the system's initial concurrence 
and the concurrence evolution of a maximally entangled state under also the one- 
sided channel. Through this channel-dependent technique one can characterize the 
entanglement dynamics under unknown channels by probing the entanglement evolution 
of a maximally entangled state alone without exploring the concrete action of the channel 
on all initial states [UJ. The entanglement factorization law [UJ has been experimentally 
verified in an amplitude decay channel [13], and in the combined channel of phase 
damping and amplitude decay [14J. The factorization law has been generalized to finite- 
dimensional bipartite system for both initially pure state[15] and mixed statepSl ITT] , 
and to determine the evolution of the G concurrence [H] of two qudits [19] . In Ref . [20J , 
the author studied the multipartite system with only one subsystem undergoing an 
arbitrary physical process and found that the entanglement evolution of the multipartite 
system can be determined uniquely by a single function of the quantum channel alone, 
irrespective to the number of qudits in the system. Therefore, in this case, the evolution 
of multipartite entanglement is in exactly the same way as that of bipartite (two qudits) 
entanglement [20] . 

In many practical situations, each particle of an entangled system is coupled with 
a local noisy environment, therefore the factorization law of entanglement evolution 
under one-sided noisy channel should be generalized to the case of multi-sided channels. 
It would be of great interest if the entanglement evolution under multi-sided noisy 
channels could be factorized as that in the one-sided channel. Unfortunately, however, 
even for the simplest case of two-qubit entanglement under two-sided noisy channels, 
the factorization law [TTJ does not hold any more except for some special forms of 
two-qubit state under special environments [21, 22J. It was shown that [2TJ for the 
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two-sided amplitude damping channels, only for the NOE (the number of excitations 
in the system is not more than one) state the factorization law is valid. For the two- 
sided pure dephasing channel case, the factorization law is valid not only for the NOE 
state but also for all the three basis states [22]. Therefore, a general factorization law 
of entanglement evolution does not hold for multi-sided channels. However, it is still 
useful to explore the factorization law for special forms of entangled states, especially 
the multipartite entanglement, under special noisy channels. This study may shed some 
light on the understanding of entanglement dynamics of multipartite system. In this 
work, focusing on the iV-qubit generalized Greenberger-Horne-Zeilinger (GGHZ) state 
[23], we shall investigate the possibility and conditions to factorize the entanglement 
evolution of multiqubit system under multi-sided noisy channels. 

In studying the dynamics of multipartite entanglement, one of the biggest obstacles 
is the lack of a computable entanglement measure. Most studies [211 ESI EHl E7J EE] on 
the entanglement dynamics of multiqubit system are based on the strategy of bipartition 
and the measure of bipartite entanglement. Another strategy one usually adopts is to 
generalize the bipartite entanglement measure directly to the multipartite case [29, 30J. 
However, the calculation for the mixed state of a multipartite system requires an 
optimization process, which is very difficult to solve exactly and can only be determined 
purely algebraically in the regime where the mixing is moderate [29]. To be able to 
assess the global entanglement of a multiqubit system, we adopt the concept of lower 
bound of concurrence (LBC) [311 [32] . Based on the derived formula of LBC evolution 
of a iV-qubit system under some typical noisy channels, we find two kinds of sufficient 
conditions under which the factorization law holds, namely, the time- dependent LBC 
can be determined by the product of initial LBC of the system and the LBC evolution of 
a maximally entangled iV-qubit state under the same multi-sided noisy channels. Then 
we consider the realistic situations under which the two kinds of factorization conditions 
can be satisfied. The dependence of entanglement robustness on the number of qubits 
and that of the noisy channels is also discussed. 

2. LBC evolution and the conditions for its factorization 

2.1. Time- dependent LBC of a N-qubit system under multi-sided noisy channels 

In this work, we adopt the LBC proposed in [32] to quantify the global entanglement of 
a multiqubit system, defined for an arbitrary iV-qubit mixed state p as 



i>l 

In Eqs. ([1]) and (J2]), k\N — k stands for the bipartitions with k qubits in one block 
and N — k ones in another. In view of the distinctions of the iV qubits, we use 




(1) 



with 
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V[k\N — k] to specify a concrete combination of k and N — k qubits in constituting 
the bipartitions k\N — k. Thus, Y2v stands for the summation over all possible concrete 
bipartitions of V[k\N - k}. In Eq.©, A^'jf ~ k] are the eigenvalues, in decreasing order, 
of the non-Hermitian matrix p(Lf i © L^ k ) p* : (Lf © L^ k ) with {L^; l x = 1, 2...ki} and 
{Lj^~ k ; l 2 = l,2...k 2 } the generators of group SO(2 fc ) and SO(2 7V_fc ) acting on the k and 
N — k qubits of a concrete bipartition V[k\N — k]. Surely, C N (p) > signifies that p 
is entangled and a separable state p always has C N (p) = 0. Yet, C N (p) = does not 
necessarily imply separability of p. The LBC of three-qubit X states was analyzed in 
Ref. [33] to demonstrate when it goes to zero. 

It is known the forms of multiqubit entanglement are diverse. In this work, the 
consideration is restricted to the GGHZ state of N qubits in the form 

1,2. ..JV ~~ a I* 1 ' ■■^ N )l,2...N + P ^2, ■•■ijv) 1 2 N , (3) 

where a,/3 G C satisfying |a| 2 + |/3| 2 = 1, ii,i 2 ...i N G {0,1} and i% = 1 ©ii,i 2 
I © i 2 ...i N = ~\- @ in with © being an addition mod 2. Consider the iV qubits are 
locally subjected to M < N independent noisy channels without mutual interactions 
between the qubits. The dynamics of the j'-th qubit is governed by a master equation 
that gives rise to a completely positive trace-preserving map (or channel) £j describing 
the evolution as pj(t) = SjPj(0), where Pj(0) and Pj(t) are, respectively, the initial and 
evolved states of the j-th qubit. Under the actions of M independent noisy channels, 
the density matrix of the initial state ([3]) of the iV-qubit system 

pjv(0) = |* ) ll2 ...jv (*o| 

= W\ 2 1*1, «2, -%)i, 2 ..jv (*i)*2, -%| + |/3| 2 |«i,«2,— *iv) 1)2 n (h,H,-^n\ 
+ a(5* \i x ,i 2 , -iN) li2 ...N (*i)*2, -«jv| + |*1, «2, •••^v) 12 _ Ar («1,«2, ■•■^] 

will evolve into a mixed state pjv(t) given simply by the composition of the M individual 
maps: 

p N (t) = £i£ 2 ...£ M p N (0). (5) 

In this work, we shall consider several paradigmatic types of noisy channels, such as 
the amplitude-damping (AD), depolarization (D) and phase-damping (PD) (dephasing) 
channels. Under the action of any one of these channels, it never appears new off- 
diagonal terms in pjv(t) but multiply the initial off-diagonal element of Pat(0) (j4]) by a 
time- dependent factor, whereas the diagonal elements of Pat(0), may give rise to new 
diagonal terms [27]. For convenience, in the representation spanned by the iV-qubit 
product states {|00...00), |00...01), |ll...ll)}, we label the 2^ diagonal elements of 
Piv(t) as a±(t), a m (t), a 2 jv-i (t),b 2 N-i (t), b m (t), bi(t) (from the top left one to 
the lower right one) and the two off-diagonal elements as d m (t) and d* m {t). Here, we 
specify the elements a m (t) and b m (t) (d m (t) and d^t)) as the time-evolution of the 
initial nonzero diagonal (off-diagonal) elements of Piv(0). The other ones a m '(t) and 
b m '{t) (with m! ^ m G {1, 2..., 2 JV ~ 1 }) are the new diagonal elements derived in the 
time evolution. Also, note that the two diagonal terms of pjv(t) corresponding to the 
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elements cii(t) and (i = 1, 2...2 Ar_1 ) have completely opposite marginals for all the 
iV individual qubits. As an example, suppose the initial density operator of a three- 
qubit system is p 3 (0) = |a| 2 |000)(000| + |/?| 2 |111)(111| + a/3*|000)(lll| + a*/3|lll)(000|, 
the evolved density operator of which is p 3 (t) = ai(t)|000)(000| + a 2 (t)|001)(001| + 
a 3 (t)|010)(010| + a 4 (t)|011)(011| + 6 4 (*)|100)(100| + 6 3 (t) 1 101 > < 10 1 1 + fe 2 (t)|110)(110| + 
6i(*)|111><111| + di(t)|000)(lll| + dj(t)|lll)(000|. In this case, the time-evolution of 
the initial nonzero diagonal (off-diagonal) elements are obviously a,i(t) and b\(t) (di(t) 
and dl(t)). The other diagonal elements a m /(t) and b m >(t) with m' = 2,3,4 are the 
new diagonal elements derived in the time evolution. The two diagonal terms, such as 
1 001) (001 1 and 1 110) (110| with the coefficients a^{t) and 62 (t) (the same for the other 
three pairs of diagonal terms with the coefficients aj(t) and b{(t), i = 1,3,4), involve 
completely opposite reduced states for all the three individual qubits. 

By virtue of Eq. (j2J), we obtain the concurrence of a concrete bipartition V[k\ N — k] 
of the evolved state /Ojv(^) as 



where d m (t) is the time-evolution of the initial off-diagonal element of the GGHZ state 
(TJJ, while a m r(t) and b m >(t) (with m! 7^ m) are the diagonal elements derived in the 
evolution. The value of m! G {1, 2...2 JV ~ 1 } is determined by the definite bipartition 
V[k\N — k]. By virtue of Eqs. (0Q) and (jHJ), the LBC of the iV-qubit system can be 
expressed as 



So far, we have derived a general formula <^ of the time- dependent LBC of a N- 
qubit system under multi-sided noisy channels. It should be pointed out that we have 
not solved the state evolution equation of the qubits' system in the sense that we do not 
know the explicit forms of the matrix elements of the evolved density operator pjv(i)- 
Actually, we are aware of the structure of the evolved density operator pjy(t) since we are 
restricted to the GGHZ initial state of the N qubits and to the special noisy channels, 
which ensures the derivation of the formula (fTl) . Inspections of the formula (JTJ) show 
that to make the factorization law hold for the LBC evolution of the iV-qubit system 
under multi-sided noisy channels, the associated bipartite concurrence C" p ^ N ~ k }(piy(t)) 
(J6]) should satisfy either one of the following two kinds of conditions. 

2.2. The first kind of condition 

An obvious situation under which the LBC ([7]) can be factorized is that the product of 
a m '(t) and b m i(t) in the formula (JHJ) is equal to zero, i.e., a m i{t)b m i{t) = Vm' 7^ m G 
{l^...^" 1 }. In this case, we have C v ^ N ~ k \p N {t)) = 2\d m {t)\ which is independent of 
the concrete bipartition V[k\N — k}. Accordingly, the LBC (JTJ) of the iV-qubit system 
is reduced to a simple form as 




(6) 





C N {p N {t)) = 2\d m (t)\, 



(8) 
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in which |d m (t)| is obviously the time-evolution of the off-diagonal element |<i m (0)| = 
| aj3 1 of the iV-qubit system. That is, the LBC of the iV-qubit system can be determined 
completely by the time-evolution of the off-diagonal element |d m (£)|. As mentioned 
above, application of any one of these noisy channels to the iV-qubit system will 
multiply the off-diagonal element of Piv(O) (j4j) by a time-dependent factor, namely, 
d m {t) = d m (fS)V{t). Therefore, the LBC (jHJ) can be further reexpressed as 

C N (p N (t)) = C N (p N (0))\V(t)\ (9) 

where C N (p N (0)) = 2\d m (0)\ = 2\aP\ is the initial LBC of the GGHZ state ©. By 
virtue of Eq. (Q, for the initially maximal entangled state of the iV-qubit system with 
C n (pn(Q)) = 1, we have C N (p^(t)) = [D(t)\, which implies [D{t) \ is the LBC evolution of 
the maximally entangled state of the iV-qubit system under the same multi-sided noisy 
channels. Therefore, the LBC evolution of the iV-qubit system under multi-sided noisy 
channels can be determined by the product of the initial LBC and the LBC evolution of 
a maximally entangled state. Actually, as we shall show in next section, the factor T>(t) 
can be determined by the parameters of the independent noisy channels that act on the 
qubits. Of course, the factorization expression Eq. (Q is conditioned on the condition 
a m >(t)b m /(t) = in (jSJ). In the next section, we shall show that this condition can be 
satisfied in AD and PD channels. 



2.3. The second kind of condition 

In some situations, the first kind of condition a m /(t)b m i(t) = Vm' ^ m e {1, 2..., 2 N ~ 1 } 
may do not hold and the LBC evolution cannot be factorized as the expression in Eq. 
09]). In this case, if \/ a m t{t)b m >(t) can be decomposed to \d m (0)\\J-' m i(t)\, the concurrence 
(EJ) of a definite bipartition V[k\N — k] will take the form 

C v ^ N - k \p N (t)) = 2\d m (0)\Q VWN - k] (t), (10) 

with Q-Pm-k]^ = max { , \V(t)\- \F m '(t)\}. By virtue of Eq. ©, LBC of the iV-qubit 
system is thus reduced to 



C N (p N (t)) =C N (p N (0)) J^zr^iE^ 2 ^ 1 ^ 1 ^)) 2 , (11) 

with C N (p N (0)) = 2|d m (0)| = 2\a/3\ being still the initial LBC of the GGHZ state 
([3]). From Eq. for C n (pn(0)) = 1 corresponding to an initially maximal 

entangled state of the iV-qubit system, we obtain its the LBC evolution as C N (pN{t)) = 
yj ' 2 n-i_ 1 X]-p(Q^^^ jV ~ fc K^)) 2 - Therefore, the LBC evolution of an arbitrary iV-qubit 
system under multi-sided noisy channels can still be determined by the product of the 
initial LBC and the LBC evolution of a maximally entangled state. In the next section, 
we shall show that the LBC evolution of the iV-qubit system under D channels can 
satisfy the above condition and exhibit the factorized form in (jlip . 
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3. The factorization LBC evolution in various noisy channels 

We have shown theoretically that under certain conditions the LBC evolution of a N- 
qubit system under multi-sided noisy channels may be factorized in such a way that 
the LBC is determined by the product of the initial LBC and the LBC evolution of a 
maximally entangled state under the same multi-sided channels. Now, an immediate 
question arises that in what realistic situations these conditions can be satisfied. In 
the following, we provide an answer to this question by analyzing the evolution of the 
iV-qubit system in the AD, D and PD channels, respectively. 



3.1. Amplitude- damping channel 

At first, we consider the AD channel which correspond to the zero-temperature 
dissipative reservoir. The factorization law of two-qubit under one-sided noisy channel 
pTj has been verified in the AD channel [T3| [H]. The action of an AD channel j for a 
qubit is described by a map £A D : p{t) = £f D p(0) with p(0) and p(t) the initial and 
evolved density matrixes of the qubit. During the time evolution a qubit decays from 
its excited state |1) to ground state |0) by emitting an excitation, with a probability 
Pj(t) = 1 — exp(— Vjt), where Tj is the decay rate of the noisy channel. The action of 
£A D on elements of the reduced density matrix of a qubit reads: 

|o> (o|-ho) <o|, 



AD 



(12) 



|0)(l|-> y /l- Pj (t) \0)(l\, 
|1> <0|->v/l-p-(f)|l> (0|, 
|l)(l|^ Pj (t)|0)(0| + (l- Pj (t))|l)(l|. 

Different forms of initial GGHZ state (J3J), such as \tp s ) = a|000) + /3 1 1 1 1} and \ip a ) = 
a|001) +/?|110) for a three-qubit system, will exhibit different dynamical evolution in the 
AD channels [28]. Therefore, we classify the initial GGHZ state of the iV-qubit system 
into symmetrical and asymmetrical cases: in the symmetrical GGHZ state (1131) all the 
iV qubits are equivalent with the same reduced states, while in the asymmetrical state 
(II 6p the N qubits can be divided to two asymmetrical blocks with opposite reduced 
states. 

The symmetrical GGHZ state can be expressed as 

l*l)l,2..JV = « l°> 0, °)l,2...7V + £ |1, 1, -l)l,2..JV ( 13 ) 

Without loss of generality, suppose the qubits 1, 2, M (M < N) are locally coupled 
to M independent and different AD channels, respectively. Under the actions of the 
M AD channels the initial density matrix p l N (0) = |^i) 12 at(^i| wm ma P on to p^it) 
which can be obtained as 

M 

P l N (t) = lit 1 - VM\H<*F |00-0)i,2...iv (11-11 + |H-l)i l2 ...iv (00-01) 

3=1 

1 

+ W ?l(*)-9Af (*) Kl-ijW, l-l)l...M,M+1...7V 1 -1| 

ii...iM=0 
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+ |a| 2 |0...0) UJV (0...0| ! (14) 

where the coefficients qj(t) = Pj{t) Wij = and qj(t) = 1 — Pj(t) Vi, = 1 for j = 
1,2,..., M. Recall that we specify a m (t) and b m (t) (d m (t) and d* m {t)) as the time evolution 
of the nonzero diagonal (off-diagonal) elements of initial density operator /Ojv(0), while 
a m i (t) and b m i (t) (with m' ^ m) are the derived diagonal elements in the evolution. From 
(JHJ) we can see that the initial off-diagonal element |d m (0)| = \a(3\ has been evolved to 
l^mCOl = \ a P\\ YljLiil ~ Pj{t)]^\i while the initial diagonal elements a m (0) = | or | 2 keep 
invariant and b m (0) = \(3\ 2 has been evolved to b m {t) = |/3| 2 (1 — pi(t))...(l — puif))- 
To decide if the LBC of p l N (t) can be factorized, we should check if the aforementioned 
condition a m ,(t)b m <(t) = 0Vm'^me {1, 2..., 2 JV ~ 1 } can be satisfied. From the second 
term of the right-hand side (RHS) of Eq. (114p . we observe that if not all the N qubits 
are subjected to the AD channels, i.e., M < N, all these diagonal terms have the same 
reduced density operators |11...1) M+1 M+2 N (11. ..1| for the qubits M + 1, M + 2...N 
that are not subjected to the noisy channels. In other words, the matrix p l N (t) (fl4l 
does not involve the derived diagonal terms whose reduced density operators for the 
qubits M + 1,M + 2,...,N are |00...0} M+1M+2 >>JV (00...0|. Recall that the elements 
a m r(t) and b m '(t) are the derived diagonal terms that should possess opposite reduced 
density operators for all the N individual qubits, therefore in the case of M < N, we 
always have a m i(t)b m i(t) = 0. The LBC of the iV-qubit system can thus be obtained by 
virtue of Eq. (EJ) as 

M 

= 2\aP\H[l- Pj (t)]t- ( 15 ) 

3=1 

Obviously, the 2\a/3\ in ffl5l) is the initial LBC of the GGHZ state (Tl3|) while the term 
n^li(l — Pj{t))* comprises the parameters of the M AD channels. 

In contrast to the case of M < N, if all the iV qubits are coupled to the AD 
channels, i.e., M = N, then all the 2 N diagonal elements of the states (IT41) are non- 
zero, which implies that the condition a m i(t)b m >(t) = cannot be satisfied. Therefore, 
the LBC of the system cannot be factorized to the form ffl5|) . Actually, it has been 
shown [27] under the action of iV independent AD channels, the iV-qubit system in the 
GGHZ state (fT3|) will suffer from the entanglement sudden death [7] when \/3\ > \a\. 
At the same time, the robustness of the iV-qubit system, i.e., the time at which the 
entanglement become arbitrarily small, decrease with the increase of the number N of 
qubits [27]. Here, for the case of M < N, we can see from the factorization expression 
( fl5|) that the LBC is independent of the system size iV and only related to the number 
M of the AD channels. In Fig. 1, we plot the evolution of LBC under M = 1, ...,5 
identical AD channels for the systems with N > M qubits. For convenience, we have 
parameterized the time dependence in terms of p = p(t) instead of t noticing that p = 
when t — and p — > 1 when t — >■ oo, i.e., p £ [0, 1] for t £ [0, oo]. From Fig.l, we can see 
that the decay rate of the system increase with the number M of the AD channels. It 
deserves noting that for M = 1, i.e., under one-sided channel, the LBC evolution of any 
system with N > 2 qubits are equivalent to the concurrence evolution of a two-qubit 
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p 



Figure 1. The LBC evolution of a system with N > M qubits that are initially in 
the symmetrical GGHZ state (fT3|) under actions of M independent AD channels. The 
parameter of the initial GGHZ state is \a\ — |/3| = 1/^/2 and p G [0, 1] for t 6 [0, oo]. 



system under also the one-sided channel, which is consistent with the result of Gour in 

Next, we investigate the asymmetrical GGHZ state in which the reduced states 
of the N qubits are not equivalent. Without loss of generality, we assume the reduced 
states of the qubits 1, 2, n with 1 < n < N are all |0), whilethe qubits n+1, n+2, N 
are all |1), therefore the second class of GGHZ state can be expressed as 

1*11)1,2...* = <> I0-..0, l...l> : ..,,,.. : v + 0...0), .,,,..,... v . (16) 

Suppose all the N qubits are coupled with their own AD channels. Then the initial 
density matrix p^(0) = l^n)^ n (*nl wm ma P on ^° Pn^) which can be obtained as 

TV 

($(t) = II^ 1 - PMk<*P* I0.-.0, 1-1)!...^+!...^ (1-1, 0...0| + h.c.) 
1 

+ \a\ 2 q n+1 (t)...q N (t) \0...0,i n+1 ...i N ) 1 n n+1 N (0...0 1 i n+1 ...i N \ 

i n+ i...i N =0 
1 

+ \(3\ 2 /;.../,..()...()),...„,,. ,... v (/,.../».()...() • (17) 

h...i„=0 

where the coefficients qj(t) = Pj(t) Vij = and qj(t) = 1 — Pj(t) Wij = 1 with 
j G {1, 2, iV}. Here, we directly check if the first kind of condition a m i{t)b m i{t) = 
VW 7^ m G {1, 2..., 2 iV-1 } can be satisfied by p l ^{t). From Eq. f JT7|) . we observe 
that there exist only two diagonal terms, i.e, |0...0, l...l) 1 nn+1 N (0...0, l...l| and 
|l...l, 0...0) 1 nn+1 N (l-.-l, 0...0|, which involve completely opposite reduced states for 
all the N individual qubits. However, these two terms are merely the nonzero diagonal 
terms of the initial density operator p l ^(0) rather than the derived terms in the 
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evolution. That is, the derived diagonal terms with elements a m t(t) and b m i(t), such 
as |0...0, l---0)i.„ nn +l,...N (0...0,1...0| and |l...l, 0...1) i...n,n+i,...jv(l---l 5 — X | , which involve 
completely opposite reduced states for all the N individual qubits will not appear 
simultaneously (for the present example only the former term can appear while its 
counterpart does not). Therefore, the condition a m i(t)b m >(t) = can be satisfied for 
the iV-qubit system in the GGHZ states ffl6|) under the action of N independent AD 
channels. By virtue of (j5J), we can get the factorized form of LBC for the state p"(i) as 

JV 

C N (p*(t)) = 2\ap\l[(l-p j (t)) 1 *. (18) 

3=1 

Obviously, 2\aft\ in (fl8|) is the initial LBC of the state (fl6|) while the term n^ =1 (l — 
Pj(t))v reflects the actions of the N AD channels. Here we have considered all the N 
qubits are subjected to the independent AD channels, i.e., M = N. For the case of 
M < N, the factorized form in (|T8|) still hold but the RHS of (fl8l) should be replaced 
by 2\a(3\ Ylj=i{l ~ Pj(X)) 5 implying the action of M AD channels. 

3.2. Depolarization channel 

The depolarizing channel j describes the situation in which a qubit remains untouched 
with probability 1 —pj(t), or is depolarized (white noise), i.e., its state is taken to the 
maximally mixed state, with probability Pj(t). The action of the map £f of the D 
channel j on elements of the reduced density matrix of a qubit reads: 

| > (o|^(i-^l)|o) (o| + a£|i><i|, 

|1> (l|->(l-Ed&)|i) (1| + ^|0) (0|, 

|o> (i\^(i- Pj (t))\n)(M, 

|1) (0H(1- Pi (t))|l)(0|. 

Without loss of generality, we suppose the qubits 1, 2, M are coupled to M (M < N) 
independent D channels, respectively. Then under the actions of the M different D 
channels the initial density matrix Ptv(0) (jlj) of the iV-qubit system will map onto pjv(t) 
which can be obtained as 

M 

l 

+ l«| 2 9l(*)-?Af(*) Nl---*M^M+l---«Ar)i...M,M+l...iv(*l-- i M 5 *M+l--iiv| 

i' 1 ...i' M =0 

1 

+ 1/3 1 2 Y] gi(t)...g M (t) i v ..i M ,i M+1 ...l N ) (z 1 ...T M ,i M+1 ...l2i)) 

H---4r=0 



£ D 



(19) 



l...M,M+l...iV 

Pj{t) 



where the coefficients qj(t) = 1 — ?4f^ when i'- = ij {il- = ij © 1), gj(t) = 1 — 

(^2 ) wnen ^ = h (i'j = h © 1) "with j = 1,2...M. Before discussing the general 
case, we first take a three-qubit system in the initial state ps(0) = |ct| 2 1000) 123 (000| + 
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|/3| 2 |111)(111| +a/3*\ 000) 123 (111 | +a*/3\ 111) 123 (000| as an example and suppose only the 
former two qubits are subjected to the D channels. By virtue of Eg. (120]) . the evolved 
density operator read 

2 

P3{t) = - Pi(*))(«0* I 00 °)l23 ( m l + <**P I m )l23 (000|) 

3=1 

+ |„|»[(i - - P^))|000>(000| + ^^|H0)(110| 

+ (i - 5^)^|oio><oio| + - ?f>)|ioo><ioo|] 

+ \m - - ^)iiii><iui + ^^iood(ooii 

+ (i - E^)^|ioi>(ioi| + E^(i - ^)|oii)(oii|]. (2i) 

From Eq. (1211) . we can see that any one of the former four diagonal terms with a 
common coefficient \a\ 2 has a counterpart in group of the latter four terms with a 
common coefficient |/3| 2 , such as |110)(110| and |001) (001 1 , so as to these two terms 
have completely opposite reduced states for all the three qubits. Therefore, the second 
factorization condition, i.e., \J a m i (t)b m >(t) can be decomposed to |rf m (0)||J- m /(t)|, is 
satisfied for the present instance and we have \/ a m > (t)b m > (t) = \a/3\\J- m i(t)\ with 
Fm>(t) = Pf^f) , (1 - £2M)I^& or 22^(1 - SM) being determined by the concrete 
value of m'. This specific example can be generalized to the general case. Actually, for 
the evolved density operator ( 120]) . we observe that if M < N the second factorization 
condition hold. In this case, we have a m i{t) = \a\ 2 qi(t)...qM(t), b m >(t) = \(3\ 2 gi(t)...gM{t), 
and y/a m >(t)b m >(t) = \ap\\jF m >{t)\ with \JF m >{t)\ = y/qi(t)...q M (t)gi{t)...g M (t). In 
addition, we note that the two diagonal elements a m r(t) and b m >(t) are equivalent so 
that \T m >(t)\ = \qi(t)...q M (t)\ = \gi(t)...g M (t)\. By virtue of Eq. (jTT]), the LBC of the 
state p®(t) can be expressed in the form 

C N {p D N (t)) = 2\a^\^^^J2( QVm ' k] ^^ ( 22 ) 

where 2\a/3\ is the initial LBC of the GGHZ state (J3J and Q^ k \ N - k \t) = 

max{0, | rij=i(l — Pj{t))\ ~ I Ylj=i ^ s f° r ^ ne case of M = A^, the second kind of 

condition does not hold any more, therefore the LBC evolution of the iV-qubit system 
cannot be factorized as the form (122]) . 

In the following, we make a brief discussion on the relations between qubits' LBC 
and the system's size N as well as the number of the noisy channels M. From (1221) . we 
conclude that in the time-evolution of the iV-qubit system, the LBC increases with the 
qubit's number iV for the fixed number M(< N) of the D channels, while decreases with 
the number M for the fixed qubit's number N. In Fig. 2 (a), we have shown the evolution 
of LBC of a four-qubit system in the GGHZ state (TT3]) under M = 1,2,3 identical D 
channels, where we can see the LBC decreases with an increase of M. Fig. 2 (b) is for 
the evolution of LBC for the systems with N = A, N = 3 and N = 2 qubits under 
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Figure 2. (a) The time-evolution of LBC for a four-qubit system in the GGHZ state 
(HH under M — 1 (solid line), M = 2 (dashed line) and M = 3 (dotted line) D 
channels, (b) The time-evolution of LBC for the system with N — 4 (solid line), 
N = 3 (dashed line) and N = 2 (dotted line) qubits under a single D channel, (c) The 
time-evolution of LBC for the system with N — 4 (solid line), N — 3 (dotted line) 
qubits under two D channels. The parameters used are \a\ = \ f3\ — l/y/2 and p £ [0, 1] 
for t € [0, oo]. 

one-sided D channel, from which we can see the robustness of LBC increases with N. 
In Fig.2(c) we plot the evolution of LBC for the system with N = 4 and N = 3 qubits 
under two identical D channels, which still shows that the robustness of LBC increases 
with N. Here, we have parameterized the time dependence in terms of p = p(t) instead 
of t noticing that p = when t = and p — > 1 when t —> oo, i.e., p e [0, 1] for t G [0, oo]. 

3.3. Dephasing channel 

The PD (or dephasing) channel represents the situation in which there is loss of quantum 
coherence with probability p(t), but without any energy exchange. The action of the 
map Ej D of the PD channel j on the density matrix of a qubit multiplies the off-diagonal 
elements by the factor (1 —p(t)) while remaining the diagonal elements invariant. Since 
the PD channels cannot lead to the population evolution of the iV-qubit system, both 
the elements a m >(t) and b m >(t) representing the derived diagonal terms in the evolution 
are zero. This implies that the first kind of condition for the factorization of the LBC of 
the iV-qubit system can always be satisfied when the iV-qubit system are coupled with 
independent PD channels. Therefore, the LBC of the iV-qubit system after evolution 
has the factorization form in Eq. (Q. 
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4. Conclusion 

In conclusion, by using LBC as a measure of multiqubit entanglement, we have explored 
the possibility to generalize the factorization law of two-qubit under one-sided channel to 
multiqubit under multi-sided channels. Instead of a general answer to this issue, we have 
considered a special form of multiqubit entanglement, i.e., the GGHZ state (j3J), and three 
typical types of noisy channels, i.e., the amplitude-damping (AD), depolarizing (D) and 
phasing damping (PD) (dephasing) channels. The explicit formulae (JTj) for the evolution 
of LBC as well as the associated bipartite concurrence ([6]) for an arbitrary bipartition are 
derived, based on which we observe that under two kinds of conditions the LBC can be 
factorized. That is, the LBC evolution of the iV-qubit system can be determined by the 
product of its initial LBC and the LBC evolution of a maximal entangled state under the 
same multi-sided noisy channels. We then study the realistic situations in which these 
two kinds of conditions can be satisfied. We have shown that (i) the iV-qubit system in 
symmetrical GGHZ state (1131) can satisfy the first kind of factorization condition when 
the number M of the independent AD channels is lesser than the number N of the 
qubits, while the iV-qubit system in asymmetrical GGHZ state (fl6l) can always satisfy 
the first kind of condition irrespective of the number M of the AD channels; (ii) in the 
D channels, the iV-qubit system in the general GGHZ state ([3]) can satisfy the second 
factorization condition when the number M of the D channels is lesser than N; (iii) 
the PD channels can always lead to a factorization of LBC of the iV-qubit system in 
GGHZ state ([3]) with a satisfaction of the first kind of condition. By virtue of the concise 
expressions of the LBC under different noisy channels, we have discussed the dependence 
of the entanglement robustness on the system size in terms of N and the number M of 
the noisy channels. The GGHZ state, though only be a special class of various types of 
multipartite entanglement, have attracted extensive researches from different aspects, 
such as the preparation [31], application [35], dynamics [24"1 1251 1261 1271 1291 128] and so 
on. Therefore, the study on the factorization law of special multiqubit entanglement 
under some typical multi-sided noisy channels can not only deepen the understanding 
of multipartite entanglement dynamics but also facilitate the corresponding calculation. 

Acknowledgments 

This work was supported by National Natural Science Foundation of China under Grant 
Nos. 10947006 and 61178012, the Specialized Research Fund for the Doctoral Program 
of Higher Education under Grant No. 20093705110001, and the Scientific Research 
Foundation of Qufu Normal University for Doctors. 

References 

[1] Nielsen M A and Chuang I L 2000 Quantum Computation and Quantum Information (Cambridge: 

Cambridge University Press) 
[2] Gisin N 1996 Phys Lett A 210 151 



Equation of motion for multiqubit entanglement in multiple independent noisy channelslA 



Zyczkowski K et al . 2001 Phys. Rev. A 65 012101 
Dodd P J and Halliwell J J 2004 Phys. Rev. A 69 052105 
Roos C F et al . 2004 Phys. Rev. Lett. 92 220402 
Carvalho A R R et al . 2007 Phys. Rev. Lett. 98 190501 
Yu T and Eberly J H 2004 Phys. Rev. Lett. 93 140404; 
Yu T and Eberly J H 2006 Phys. Rev. Lett. 97 140403; 
Eberly J H and Yu T 2007 Science 316 555 
Santos M F et al . 2006 Phys. Rev. A 73 040305 

Horodecki M, Shor P W and Ruskai M B 2003 Rev. Math. Phys. 15 629 
Bai Y K et al . 2009 Phys. Rev. A 80 044301 

Konrad T, DeMelo F, Tiersch M, Kasztclan C, Aragao A and Buchleitner A 2008 Nat. Phys. 4 99 
Wootters W K 1998 Phys. Rev. Lett. 80 2245 

Farias O J, Latune C L, Walborn S P, Davidovich L and Ribeiro P H S 2009 Science 324 1414 
Xu J S, Li C F, Xu X Y, Shi C H, Zou X B and Guo G C 2009 Phys. Rev. Lett. 103 240502 
Li Z G, Fci S F, Wang Z D and Liu W M 2009 Phys.Rev.A 79 024303 
Liu Z and Fan H 2009 Phys. Rev. A 79 032306 
Yu C S, Yi X X and Song H S 2008 Phys. Rev. A 78 062330 
Gour G 2005 Phys. Rev. A 71 012318; 
Gour G 2005 Phys. Rev. A 72 042318 

Tiersch M, dcMelo F and Buchleitner A 2008 Phys. Rev. Lett. 101 170502 
Gour G 2010 Phys. Rev. Lett. 105 190504 
Li J G, Zou J and Shao B 2010 Phys. Rev. A 82 042318 
Li J G, Zou J and Shao B 2011 Phys. letts. A 375 2300 

Greenberger D M, Home M A and Zeilinger A 1989 in Bell's Theorem, Quantum Theory and 

Conceptions of the Universe (Dordrecht: Luwer) 
Simon C and Kempe J 2002 Phys. Rev. A 65 052327 
Diir W and Briegcl H J 2004 Phys. Rev. Lett. 92 180403 
Hcin M, Diir W and Briegcl H J 2005 Phys. Rev. A 71 032350 

Aolita L, Chaves R, Cavalcanti D, Acin A and Davidovich L 2008 Phys. Rev. Lett. 100 080501 
Man Z X, Xia Y J and An N B 2008 Phys. Rev. A 78 064301 
Carvalho A R R, Mintert F and Buchleitner A 2004 Phys. Rev. Lett. 93 230501 
Minter F, Kus M and Buchleitner A 2005 Phys. Rev. Lett. 95 260502 
Mintert F, Kus M and Buchleitner A 2004 Phys. Rev. Lett. 92 167902; 
Chen K, Albevcrio S and Fei S M 2005 Phys. Rev. Lett. 95 210501; 
Breuer H P 2006 Phys. Rev. Lett. 97 080501; 
Ou Y C, Fan H and Fci S M 2008 Phys. Rev. A 78 012311; 
Aolita L, Buchleitner A and Mintert F 2008 Phys. Rev. A 78 022308; 
Schmid C et al . 2008 Phys. Rev. Lett. 101 260505 
[32] Li M, Fei S M and Wan Z X 2009 J. Phys. A 42 145303 
[33] Wcinstcin Y S 2010 Phys. Rev. A 82 032326 
[34] Leibfried D et al 2005 Nature (London) 438 639; 

Lu C Y et al 2007 Nat. Phys. 3 91 
[35] Hillery M, Buzek V and Berthiaumc A 1999 Phys. Rev. A 59 1829; 

Man Z X, Xia J Y and An N B 2006 J. Phys. B: At. Mol. Opt. Phys. 39 3855; 
Man Z X, Xia J Y and An N B 2007 J. Phys. B: At. Mol. Opt. Phys. 40 1767 



